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Abstract. The Asakura-Oosawa-Vrij (AOV) model of colloid-polymer mixtures 
idealizes nonadsorbing polymers as effective spheres that are fixed in size and 
impenetrable to hard particles. Real polymer coils, however, are intrinsically 
polydisperse in size (radius of gyration) and may be penetrated by smaller particles. 
Crowding by nanoparticles can affect the size distribution of polymer coils, thereby 
modifying effective depletion interactions and thermodynamic stability. To analyse 
the influence of crowding on polymer conformations and demixing phase behaviour, we 
adapt the AOV model to mixtures of nanoparticles and ideal, penetrable polymer coils 
that can vary in size. We perform Gibbs ensemble Monte Carlo simulations, including 
trial nanoparticle-polymer overlaps and variations in radius of gyration. Results are 
compared with predictions of free- volume theory. Simulation and theory consistently 
predict that ideal polymers are compressed by nanoparticles and that compressibility 
and penetrability stabilise nanoparticle-polymer mixtures. 
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1. Introduction 

Colloid-polymer mixtures are among the most actively studied materials because of 
the windows they open onto the rich physical behaviour of both soft (macromolecular) 
and hard (atomically ordered) condensed matter. The close analogy between colloidal 
particles and atoms, combined with real-time imaging of slowly moving particles, have 
deepened our fundamental understanding of phase transitions, from freezing/melting 
to the glass transition, gelation, and demixing pQ EJ [3]. With steady advances in 
nanoparticle synthesis and characterization, much recent interest has turned to the 
physical properties of nanoparticle-polymer composites jl]. 

Adding nonadsorbing polymers to a stable suspension of colloids or nanoparticles 
can induce aggregation and demixing through a depletion mechanism first explained 
over a half-century ago by Asakura and Oosawa [5]. Depletion of polymers from the 
space between two particles creates an unbalanced osmotic pressure that drives the 
particles together. Equivalently, the larger volume available to polymers amidst particles 
whose excluded-volume shells overlap increases the polymer entropy and manifests as 
an effective interparticle attraction. Tuning the range and strength of the attraction by 
varying the polymer size and concentration directly influences phase stability. 

The most widely studied model of mixtures of particles and nonadsorbing polymers 
is the Asakura- Oosawa- Vrij (AOV) model [HI E], which treats the particles as hard 
spheres and the polymers as mutually non-interacting (ideal) effective spheres of fixed 
radius, having hard interactions with the particles. Despite its simplicity, this important 
reference model succeeds in qualitatively explaining demixing, freezing, and other 
phenomena observed in real particle-polymer mixtures. Bulk thermodynamic phase 
behaviour of the AOV model has been intensively explored using thermodynamic 
perturbation theory [7J, free- volume theory [8], and a variety of computational methods. 
Most simulation studies have been based on Monte Carlo algorithms implemented in 
either the Gibbs ensemble [91 [TOj [HJ QEl [13] or the grand canonical ensemble [HI [151 EE] > 
both of which circumvent complications associated with phase interfaces. Quantitative 
discrepancies between predictions and experiments have motivated enhancements of the 
AOV model to include more realism in the polymer properties, including polymer- 
polymer interactions [TF1 HHJ [T9J, [201 ED E21 E21 EI], polydispersity in molecular 
weight [251 ESI [27], and conformational freedom of polymers on a lattice [281 |29]. Recent 
theoretical and simulation studies of particle-polymer mixtures also have explored 
interfacial properties [301 EH E2], demixing in confinement [331 EH ESI EH E3 EH], and 
dynamical properties, such as diffusion and response to shear [23l 124] . 

A polymer coil has a size that is well characterized by its radius of gyration [39] . 
Scattering experiments (using, e.g., light or neutrons) probe the mean coil size 
averaged over an ensemble of conformations. Even for a hypothetical solution having 
uniform molecular weight (i.e., chain length for linear polymers), the (pre-averaged) 
radius of gyration exhibits a broad distribution due to statistical fluctuations in coil 
conformations [101 HTJ H2]. This intrinsic polydispersity can influence the phase 
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behaviour of particle-polymer mixtures. 

To accommodate intrinsic polymer polydispersity in particle-polymer mixtures, 
Denton and Schmidt [13] recently proposed a modified AOV model in which the polymer 
coils have a single internal degree of freedom, namely size. Taking as input the radius 
of gyration distribution of an ideal polymer coil, they developed and applied a classical 
density-functional theory, which reduces to a free-volume theory for uniform fluids. In 
the "colloid" limit, in which the polymer coils are impenetrable to the larger particles, 
the theory predicts compression of polymer with increasing particle concentration, and 
a resultant stabilisation of the mixture. 

In the "protein" (or nanoparticle) limit, the particles are small enough to 
penetrate the polymer coils. Recent experiments [HI 05] H6] H7J HHJ 09] [50], computer 
simulations [29] [51] [52] [53], and theories j5T] [54] [55] [56] have begun to explore the 
behaviour of such asymmetric mixtures. These studies raise prospects of modifying 
protein solutions by adding polymer and tuning properties of polymer-metal (or 
polymer-semiconductor) nanocomposites by adding nanoparticles to a polymer matrix. 
Depending on solvent quality and other sample conditions, experiments indicate 
compression (shrinking) [57] EH], expansion (swelling) [59], or little change in size [60] of 
polymers in response to nanoparticles. Monte Carlo simulations of bead-spring polymers 
in the presence of nanoparticles [52] [53] indicate that nanoparticles can penetrate much 
larger polymer coils. Bulk demixing behaviour has been explored via density-functional 
theory within a relatively simple model in which the nanoparticles can penetrate the 
polymers (effective spheres of fixed size) after surmounting an energy barrier [5l] . 

The main purpose of this paper is to study, by means of simulation and theory, 
the influence of nanoparticles on the conformations of nonadsorbing polymers and the 
resulting phase stability of nanoparticle-polymer mixtures. Section [2] first defines a 
simple extension of the AOV model that incorporates nanoparticle-polymer overlap 
and intrinsic polydispersity in polymer size. Section [3] next outlines our methods: 
Gibbs ensemble Monte Carlo simulation and a mean-field free-volume theory (with 
details consigned to appendices). Section 0] presents simulation results and theoretical 
predictions for polymer size distributions and demixing phase diagrams. Finally, Sec. [5] 
summarizes and concludes. 

2. Model 

We consider a mixture of impenetrable nanoparticles and nonadsorbing polymers, 
dispersed in a solvent, in osmotic equilibrium with a reservoir of pure polymer solution. 
Exchange of polymers between the system and reservoir (e.g., via a semi-permeable 
membrane) maintains constant polymer chemical potential. The thermodynamic state 
of the system is characterized by the temperature T, the nanoparticle number density 
p n , and the polymer number density in the reservoir p r p . Equality of polymer chemical 
potentials in the system and reservoir determines the polymer density in the system p p . 
To model this system, we extend the AOV model to the protein limit in the manner 
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proposed by Schmidt and Fuchs [51] by allowing nanoparticles to penetrate polymers 
and attributing to each overlapping nanoparticle-polymer pair an energy cost e. With 
the fixed nanoparticle radius denoted by R n (diameter a n = 2R n ) and the instantaneous 
radius of a polymer by R p , the nanoparticle-nanoparticle and nanoparticle-polymer 
interactions are specified by pair potentials: 

r < cr„ , 



oo, 
0, 

0. 



r>o n 
r < R n + R 



v ■ 



(1) 



(2) 



r > Rn + Rp , 

where r is the centre-to-centre distance. Assuming ideal polymers, the polymer-polymer 
interaction vanishes, i.e., v pp (r) = for all r, which strictly applies only to theta 
solvents [39], wherein the polymer second virial coefficient vanishes. 

We further extend the AOV model to describe polymer coils whose size distribution 
is intrinsically polydisperse [43] . As noted above, even an idealized solution of polymers 
with uniform chain length has a broad distribution of radius of gyration. In fact, 
the radius of gyration of an ideal, freely-jointed chain follows the exact probability 
distribution EH 1621 



Pr(Rp) 



E 



(2k + 1) 



■(4A; + 3) 7/2 exp(-t 



(3) 



where 



R r g = Jj o °°dR p RlP r (R p 



(4) 

is the root-mean-square radius of gyration of polymers in the reservoir, t — (Rp/ -^o) > 
t k = (4k + 3) 2 /(8t), and K n are modified Bessel functions of the second kind. 

The extended AOV model is fully specified by the polymer size distribution 
[equation and two parameters: the penetration energy e and the ratio of the rms 
radius of gyration of polymers in the reservoir to the nanoparticle radius, q r = R r g /R n , 
which is an experimentally accessible property. In the colloid limit (q r < 1), we assume 
impenetrable polymers by taking infinite penetration energy (e — > oo) in equation (j2J). 
In the protein limit (q r ^> 1), in which the particles can penetrate the polymers, we 
follow Schmidt and Fuchs [54] and take for the penetration energy an approximation 
from polymer field theory [63[ El] for the average excess free energy cost of inserting a 
hard sphere into an ideal coil: 

H = ~ , (5) 

q 

where = l/k^T and q = R p /R n for a polymer of radius R p . The loss of conformational 
entropy of a polymer that harbours a nanoparticle is thus modelled by an energy 
penalty that decreases as the polymer swells. (More refined models replace the step- 
function profile of equation (jSJ) by a continuous function of separation [651 66J.) Next 
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we explore the combined influences of penetration and compressibility on polymer size 
and demixing. 

3. Methods 

3.1. Gibbs Ensemble Monte Carlo Simulation 

The Gibbs ensemble Monte Carlo (GEMC) method M EH H21 E3] provides a 
computationally efficient means of computing phase coexistence curves for model fluids. 
Allowing each phase to occupy a separate box avoids any need to simulate interfaces. 
This method is known to be of limited accuracy near critical points [67J [68j [691 HO], 
where grand canonical ensemble methods [HI [151 EH] can yield higher resolution. For our 
purposes of exploring nanoparticle-polymer demixing, however, GEMC proves practical. 

We work in the semigrand ensemble, where temperature, total number of 
nanoparticles, and total volume, are fixed, while exchange of polymers with a reservoir 
fixes the polymer chemical potential. The conventional GEMC algorithm for our system 
involves four types of trial move: (1) displacements of nanoparticles and polymers within 
each box to ensure thermal equilibrium; (2) exchanges of volume between boxes to 
ensure mechanical equilibrium, characterized by equality of pressures; (3) transfers of 
nanoparticles and polymers between boxes to ensure chemical equilibrium, characterized 
by equality of chemical potentials for each species; and (4) transfers of polymers between 
each box and the reservoir. The acceptance probabilities for these standard trial moves 
are given in Appendix A. 

To explore the influence of nanoparticles on the polymer size distribution, we 
implement an additional trial move: variation in polymer radius of gyration R p . This 
move allows the polymer size distribution to adjust to the presence of particles. A 
polymer coil has radius of gyration R p with probability P(R P ) oc P r (R p ) exp(— (3U), 
where P r (R p ) is the polymer size distribution of equation ([3]) and U is the potential 
energy due to any penetrating particles. A trial change in a polymer's radius of gyration, 
from its old value R° to a new value R p , with an attendant change in potential energy 
AU, is then accepted with probability (see Appendix A) 



Through trial expansions and contractions, the polymers achieve their equilibrium 
size distribution. In a dilute suspension of nanoparticles, the size distribution approaches 
that of the polymer reservoir P r (R p ). With increasing nanoparticle concentration, 
however, crowding and penetration influence the polymer size distribution. Exploring 
the shift in polymer size distribution induced by nanoparticles, and the resulting effect 
on demixing, is the main goal of our simulation study. 




(6) 
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3.2. Simulation Details 

We simulated mixtures of nanoparticles and polymers confined to a cubic box (or 
boxes) with periodic boundary conditions. In the protein limit, we performed GEMC 
simulations in the semigrand ensemble, with fixed total number of nanoparticles N n , 
total volume V, and polymer reservoir density p p , for reservoir size ratio q r = 3. The 
volumes of the two boxes, initially set equal, were determined by N n and the nanoparticle 
volume fraction r] n = (4ir/3)p n R 3 1 . In the colloid limit, we performed canonical ensemble 
simulations (one box), with fixed N n , polymer number N p , and volume, for q r = 1 . 

The initial configurations were generated by randomly placing particles and 
polymers (initially monodisperse) on fee lattice sites. At polymer concentrations 
exceeding full lattice occupancy, some sites were multiply occupied by polymers. The 
initial total number of polymers was chosen to yield a desired effective polymer volume 
fraction, defined as r\ p = (47i/3)p p (R r g ) 3 - - an experimentally controllable quantity. 
From these initial states, the simulations proceeded via the various trial moves outlined 
in Sec. 13.11 Tolerances were adjusted to yield practical acceptance ratios for each move. 

Several diagnostic quantities were calculated to confirm equilibration. The pressure 
p in each box was computed from a simple adaptation of the virial expression for a binary 
hard- sphere fluid mixture [7T] : 

^ = 1 + \np (x 2 n a 3 n g nn (a n ) + 2(1 - e~^)x n £ x p (R n + R p fg np {R n + Rp) ) , (7) 

P 6 \ R p I 

where in the given box p is the average total number density, x n and x p are the average 
concentrations of nanoparticles and polymers, g nn {a n ) and g np (R n + R p ) are the contact 
values of the nanoparticle-nanoparticle and nanoparticle-polymer radial distribution 
functions, the summation runs over all polymer radii, and angular brackets represent an 
ensemble average over all configurations of nanoparticles and polymers. In practice, the 
contact values of the radial distribution functions were computed by assigning polymer 
radii to bins of width (R n + R r g )/50, accumulating particles in radial bins of width 5% 
of the radius of the central nanoparticle, and extrapolating to contact. 

The chemical potential of the nanoparticles was computed using Widom's test 
particle insertion method [72], applied to trial nanoparticle transfers: 



0fln = _ln (i^: exp( "^ AC/) ) ' (8) 

where AU is the change in potential energy of the box into which a nanoparticle 
is inserted and ( ) represents an ensemble average over all configurations. In the 
semigrand ensemble, the chemical potentials of polymers of all sizes are imposed by 
the reservoir. As a consistency check, the mean polymer chemical potential (averaged 
over polymer sizes) also was computed via Widom's insertion method applied to trial 
polymer transfers: 

0Mp = - ln {Z^l exp(-/3Al7)\ , (9) 
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where AU now is the change in potential energy upon insertion of a polymer. Note 
that the chemical potentials are defined relative to an arbitrary reference potential, 
determined by the choice of volume units (here a 3 ). 

For comparison with free- volume theory (Sec. I3.3[) . we also computed the effective 
fraction of the total volume (averaged over polymer sizes) that was accessible to 
the polymers, i.e., not excluded by the nanoparticles. At equilibrium, the chemical 
potentials of polymer in the system and reservoir must be equal. Thus, equating \i v 
[from equation (jHJ)] to the mean chemical potential of (ideal) polymers in the reservoir, 

Pfi* = In (p r p a 3 n ) = In (p p al/a eS ) , (10) 

the effective polymer free-volume fraction amidst nanoparticles can be expressed as 

a e ff = M^^exp(-/3A[/)\ . (11) 

The effective polymer free- volume fraction for a given box was determined by applying 
equation (ITTj) during the GEMC simulation and averaging over all trial polymer transfers 
into the box. 



3.3. Free-Volume Theory 

To guide the choice of parameters in our simulations, we generalize the free-volume 
theory of Lekkerkerker et al [8] to the case of compressible, penetrable polymers. As 
shown in Appendix B, the Helmholtz free energy density / of a mixture of nanoparticles 
of number density p n (volume fraction r] n ) and polymers in osmotic equilibrium with a 
polymer reservoir of density p r p can be expressed (in ksT units) as 

f(Vn, Pi) = Pn [in (pnO-l) - l] + (j)hs(Vn) + P^ff^n) [in (p£<) - l] , (12) 

where 

roo 

a e &(Vn) = / dR p P r (R p )a(R p ;r] n ) (13) 
Jo 

is the effective free-volume fraction of the polymers amidst nanoparticles of volume 
fraction 7] n [cf. equation fill}) ], expressed as a weighted average of the free- volume fraction 
a(R p ; 7] n ) of polymers of radius R p . The first two terms on the right side of equation ffl2l 
are the ideal-gas and excess free energy densities of the hard-sphere nanoparticles; the 
third term is the total free energy density of the polymers. Note that in this mean-field 
approximation, which neglects nanoparticle-polymer correlations, a(R p ; r] n ) and a e s(rj n ) 
are assumed to be independent of the polymer density. The corresponding polymer size 
distribution in the mixture is given by (see Appendix B) 

P(R p , Vn ) = ^f^P r (R p ). (14) 
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Figure 1. Phase diagram of a nanoparticle-polymer mixture in osmotic equilibrium 
with a reservoir of polymers whose mean radius of gyration is three times the 
nanoparticle radius (q r = 3). Symbols represent simulation data and curves predictions 
of free-volume theory for nanoparticle-rich and -poor binodals (open circles are 
predicted critical points). Results are shown for three polymer models: AOV model 
with impenetrable, incompressible polymer (black circles, solid curve); penetrable, 
incompressible polymer (blue squares, dashed curve); and penetrable, compressible 
polymer (red triangles, dotted curve). Inset: expanded view of nanoparticle-rich phase. 

4. Results and Discussion 

We initialized the GEMC simulations with N n = 400 nanoparticles and N p = 1000 
polymers in a volume yielding a nanoparticle volume fraction of r] n = 0.03. After 
equilibrating for 150,000 MC cycles, statistics were accumulated and averaged over the 
next 150,000 cycles. Equilibrium was diagnosed via equality of pressures and chemical 
potentials in coexisting phases. Several longer runs for larger systems confirmed that 
equilibrium was attained and that finite-size effects were negligible. The main simulation 
results reported below are coexistence densities of nanoparticle and polymer species and 
probability distributions (histograms) of polymer radius. 

Figure [1] shows the fluid-fluid demixing phase diagram for reservoir polymer-to- 
nanoparticle size ratio q r = 3 (nanoparticle limit). In this representation, the polymer 
volume fraction in the system (rather than in the reservoir) is plotted on the vertical axis, 
facilitating comparison with experiment. Corresponding points on the nanoparticle-rich 
and -poor branches of the binodal represent coexisting "liquid" and "vapour" phases. 
The demixing and freezing transitions being well separated for this size ratio, we need 
not consider the liquid-solid phase boundary. Results are shown for the AOV model 
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Figure 2. Probability distribution of polymer radius of gyration in a mixture of 
nanoparticles and penetrable, compressible polymers in osmotic equilibrium with a 
polymer reservoir. The mean radius of polymer in the reservoir is three times the 
nanoparticle radius (q r — 3). Polymer size distributions are shown along the liquid 
and vapour binodals for polymer chemical potentials —j3fj, p — 0.9, 1.3, 1.5, and 1.7 (left 
to right), corresponding to mean nanoparticle volume fractions r\ n = 0.120, 0.0809, 
0.0572, and 0.0312. The symbols are simulation data, solid curves are predictions 
of free-volume theory, and the dashed curve is the reservoir size distribution (pure 
polymer solution). Inset: Ratio of mean polymer radius of gyration in the system 

1/2 

R g = (Rp) [equation (TT5l) ] to radius of gyration in the reservoir R r g (see Table 1). 

and its two penetrable-polymer extensions (incompressible and compressible polymers). 
Numerical data from our simulations are plotted together with predictions from our free- 
volume theory calculations, derived from a coexistence analysis that equates pressures 
and chemical potentials in each phase [using the free energy of equation (I12j) ]. 

The simulations indicate that polymer penetrability and compressibility both 
promote stability against demixing, raising the binodal relative to that of the AOV 
model. The free-volume theory predicts the same qualitative trend, but somewhat 
underestimates stability of the mixture — a result of neglecting nanoparticle-polymer 
correlations. The theory also predicts an increase in critical polymer concentration with 
added polymer freedom. Although this prediction is impractical to test with the Gibbs 
ensemble method, it could be tested by grand canonical ensemble methods. 

Polymer radius of gyration distributions (normalized histograms) are shown in 
Fig. |2] over a range of polymer chemical potentials along the liquid-vapour binodal. 
Our simulation data (smoothed by a running average) are compared with theoretical 
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Table 1. Polymer compression ratio R g (rj n ) / BL for reservoir polymer-to- nanoparticle 
size ratio q r = R r g / R n = 3 (see inset to Fig. [2} . Results of Monte Carlo simulations are 
compared with predictions of free- volume theory over a range of nanoparticle volume 
fractions rj n and polymer volume fractions r\ v along the demixing binodal. 



Vn 




Simulation 


Theory 


0.002 


3.21 


0.987 ± 0.001 


0.991 


0.009 


2.17 


0.968 ± 0.005 


0.963 


0.013 


1.93 


0.959 ± 0.011 


0.948 


0.042 


0.731 


0.882 ± 0.011 


0.859 


0.049 


0.587 


0.863 ± 0.007 


0.841 


0.057 


0.459 


0.842 ± 0.005 


0.822 


0.069 


0.317 


0.811 ± 0.004 


0.796 


0.081 


0.227 


0.786 ± 0.002 


0.773 


0.091 


0.173 


0.768 ± 0.001 


0.754 


0.100 


0.133 


0.749 ± 0.001 


0.739 


0.110 


0.106 


0.735 ± 0.001 


0.723 


0.120 


0.082 


0.720 ± 0.002 


0.707 



predictions from equation f|T4|) . The size distribution is seen to shift towards smaller 
radii of gyration with increasing nanoparticle concentration, reflecting compression of 
polymers due to crowding and penetration by nanoparticles. Endowing each polymer 
with the internal freedom to change its radius of gyration allows the polymers to shrink 
(compress) to avoid penetration by nanoparticles and contributes to stabilising the 
mixture. Polymer compression is further quantified by the mean radius of gyration, 

R 9 (Vn) = (Rl) 1/2 = URpY, P(R^Vn)(R^ ' , (15) 

defined as an average over conformations, i.e., sum over histogram bins, where is 
the polymer radius at bin i, and AR p is the bin width. The corresponding polymer 
compression ratios R g /R r g are plotted in the inset of Fig. [2] and listed in Tabled] 

To further explore polymer compression, we performed canonical ensemble 
(constant- NVT) simulations in the colloid limit (q r = 1) at fixed polymer volume 
fraction. Polymer size distributions and compression ratios are shown over a range 
of particle volume fractions in Fig. El For dilute suspensions, crowding by particles 
compresses the polymers, as for q r — 3. These results support the trends in polymer size 
distribution predicted previously (13] for the same compressible polymer model. With 
increasing particle concentration, however, the polymers cease shrinking and begin to 
expand. This size reversal can be understood by noting that the trajectory of the state 
point now crosses the binodal from the stable (mixed) region into the unstable (demixed) 
region of the phase diagram. Snapshots of the system in the unstable region reveal the 
presence of large polymer clusters. Within such a cluster, polymers are shielded from 
particles and thus are free to adopt a size distribution closer to that in the reservoir. 
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R / R r 

p g 

Figure 3. Probability distribution of polymer radius of gyration R p in a mixture of 
particles and impenetrable, but compressible, polymer at fixed polymer volume fraction 
7] p = 0.1. The mean radius of polymer in the reservoir is equal to the particle radius 
(q r = 1). Simulation data (symbols) and predictions of free- volume theory (curves) are 
shown for increasing particle volume fraction r) n = 0, 0.1, 0.2, and 0.3 (right to left). 

Such correlation-driven behaviour clearly is not captured by the mean-field theory. 

The influence of polymer compressibility on phase behaviour can be interpreted also 
in terms of effective interactions between nanoparticles induced by polymer depletion. 
Note that compression of polymers by crowding shortens the range of depletion-induced 
attraction between nanoparticles to an extent that depends on particle concentration. 
Since polymers in the nanoparticle-rich phase are more compressed, on average, than 
those in the nanoparticle-poor phase, weakening of the effective attraction tends to 
stabilise the nanoparticle-rich phase against demixing. 

Figure H] finally presents results for the polymer free-volume fraction a, computed 
from equations (11 ip and (fl3|) . along the binodal of Fig. [IJ Polymer compression clearly 
increases the free- volume fraction accessible to the polymers, consistent with the shift of 
the nanoparticle-rich binodal towards higher polymer volume fractions. The free-volume 
theory captures the qualitative trend, but underestimates a at higher nanoparticle 
concentrations, again because of neglect of nanoparticle-polymer correlations. 

5. Summary and Conclusions 

In summary, we have performed Monte Carlo simulations of model mixtures of 
monodisperse hard-sphere nanoparticles and nonadsorbing, intrinsically polydisperse 
polymers in a theta solvent. The polymers are modeled as effective spheres, each 
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Figure 4. Effective polymer free- volume fraction a e g in a mixture of nanoparticles 
and penetrable, compressible polymer in osmotic equilibrium with a polymer reservoir 
for reservoir polymer-to-nanoparticle size ratio q r = 3. The symbols are simulation 
data [equation and the curve is the prediction of free-volume theory [Eq. (|13[) ] 

along the liquid and vapour (demixing) binodals. 

with a single internal degree of freedom (average size), which allows for fluctuations 
in radius of gyration. In addition to conventional Monte Carlo moves, the polymers 
undergo trial size variations, with an acceptance probability dependent on nanoparticle 
concentration. Within this coarse-grained model, we have investigated the effect of 
nanoparticle crowding on polymer size and the influence of polymer compression and 
penetrability on the demixing behaviour of nanoparticle-polymer mixtures. 

For particles and polymers of comparable size, our simulations confirm the trends 
previously predicted by density- functional theory [33], namely that particles compress 
the polymer coils and that mixing stability is enhanced by polymer compression. In 
equilibrium fluid-fluid coexistence, the polymers fractionate into compressed coils in 
the particle-rich phase and expanded coils in the particle-poor phase. In the protein 
limit, in which nanoparticles face an energy barrier for penetrating the larger polymers, 
our simulations also show significant polymer compression and enhanced stability, in 
qualitative agreement with predictions of free-volume theory. Incorporating into the 
model excluded-volume interactions between polymer segments may help to shed light 
on the conformations of crowded polymers in good solvents. 

While our model includes intrinsic polydispersity of polymers, it neglects 
polydispersity in chain length (molecular weight), which can be significant in real 
polymer systems [25j EEl E7]. Fasolo and Sollich [27] have applied free- volume theory 
to mixtures of particles and polydisperse polymers, finding that polydispersity in chain 
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length promotes demixing - an effect just opposite to that of polymer compressibility. 
The relative importance of these two competing effects could be assessed by combining, 
in a single model, both radius of gyration and chain length polydispersities. 

Finally, we emphasize that the present study demonstrates a general and practical 
approach to modelling soft matter systems, in which the microscopic complexity 
of the constituent macromolecules is considerably reduced through a coarse-grained 
approximation, while physically relevant features are retained and subsumed into a small 
number of internal degrees of freedom. Conceptual insight provided by this mesoscopic 
approach may help to guide the design of experiments and the choice of parameters in 
studies of more detailed molecular models [291 E21 E3]. A similar analysis of polymer 
shape variations induced by particle crowding, and the influence on demixing behaviour, 
will be the subject of future work. 
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Appendix A. Acceptance Probabilities of Monte Carlo Trial Moves 

The acceptance probabilities for the GEMC trial moves follow directly from the 
condition of detailed balance [73], according to which the average rate of transition from 
an old state (o) to a new state (n) must equal, in equilibrium, the average transition 
rate from state n to state o. Defining P{6) and P(n) as the probabilities of finding the 
system in the states o and n, respectively, and 7r(o — > n) as the transition probability 
between the states, detailed balance requires that 



Assuming that trial moves o — > n and n — > o are attempted at equal rates, equation 
flA.lj) implies acceptance probabilities in the ratio 

accp ->• n) = P{n) 
acc(n -»■ o) P(o) 

In the Metropolis algorithm, trial moves are accepted with probability [731 E] 



Trial displacements of randomly chosen particles are thus accepted with probability 




(A.l) 




(A.3) 



P disp = min{l, exp(-/3AC/)} , 



(A.4) 
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where AU = U n — U a is the change in energy between old and new configurations. A 
trial move leading to overlap of nanoparticles yields infinite AU and so is automatically 
rejected. Each nanoparticle-polymer overlap increases U by e; otherwise AU = 0. 

In the semigrand Gibbs ensemble, box % holds Ni particles in a volume Vi (i = 1, 2). 
A trial exchange of volume AV, achieved by uniformly rescaling all particle coordinates, 
such that Vi — > V\ + AV and V 2 — > V 2 — AV, is accepted with probability 

(Y^f (liz^f exp( _, At/) } . (A . 5) 

In practice, it proves more efficient to make trial moves in \n(Vi/V 2 ), for which the 
acceptance probability is [731 EI] 

f /Vi + AV\ Nl+1 fVo - AV\ N2+1 1 
P vol = min|l, [ y ) ( 2 y ) exp(-/3AC/)| . (A.6) 

Denoting by ATji and the respective particle numbers of species j (j = n,p) 
in the two boxes, transfer of a particle of species j from box 1 to box 2, resulting in 
Nji — > Aji — 1 and Nj 2 — >■ A"j 2 + 1, is accepted with probability [H] 

n rans = min jl, ^M^exp(-/3Af/)| . (A.7) 

In the colloid limit (g r < 1), because of the difficulty of inserting a large particle into a 
box without overlaps, only direct transfers of polymers between the two boxes were 
attempted. Transfers of particles were achieved indirectly by exchanging identities 
of particles and polymers p2J [13]. In an identity-exchange move, a randomly chosen 
polymer in a randomly chosen box is changed to a particle; in the other box, a randomly 
chosen particle is changed to a polymer. Changing a polymer to a particle in box 1 (and 
a particle to a polymer in box 2) is accepted with probability 

v - = miD {^ jK^wt^) eM ^ AUl+AU2>] } • (A ' 8) 

where N nl (N p i) and N n2 (N p2 ) are the respective numbers of particles (polymers) in 
the two boxes and AU\ and AU 2 are the resulting changes in energy in boxes 1 and 2 
- incremented by e for each overlap of a particle and a polymer. 

Trial transfers of polymer between the system and the reservoir are attempted with 
equal frequencies in either direction. Transfer of a polymer from the reservoir to the 
system is executed as follows. A polymer with a radius chosen randomly from the 
reservoir size distribution is placed at a random position in a randomly chosen box. For 
a transfer from the reservoir to box 1, the acceptance probability is 

J\^Ti eXp( ~ /9At/l) } ' (A ' 9) 
where AU\ is the change in energy of box 1. Transfer of a polymer from the system to 
the reservoir proceeds by removing a randomly chosen polymer from a randomly chosen 
box, with acceptance probability 

P res = min|l, NpX ~ l exv{-PAU x )\ . (A.10) 
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Finally, from equation (I A . 3 j) . a trial change in polymer size from radius R° to R™ 
is accepted with a probability given by equation fl6]). The same result follows also from 
the general relation V(o — > n) — min {1, exp(— (3AF)}, where AF is the change in 
Helmholtz free energy, if the conformational entropy of a polymer in the system is taken 
to be the same as in the reservoir, i.e., — k^Thi P r (R p ). 



(f>hs{Vn) = Pn— ZTTT~ ■ ( B ' 2 ) 



Appendix B. Free- Volume Theory of Nanoparticle-Polymer Mixtures 

Following [HI H3] , the Helmholtz free energy density / of a mixture of nanoparticles of 
number density p n (volume fraction r] n ) and polymers in osmotic equilibrium with a 
polymer reservoir of density p r p can be expressed (in ksT units) as 

f(Vn, Pp) = Pn [in (pnO-l) ~ l] + fihsiVn) + fp(Vn, Ppl , (B.l) 

to within an arbitrary constant. The excess free energy density of the hard-sphere 
nanoparticles is accurately approximated by the Carnahan-Starling expression |71j : 

r) n (4 - 3r ] n ) 

i 1 ~ Vn) 

Progress in approximating the polymer free energy density f p is facilitated by defining 
polymer density distributions p p (R p ;ri n ) and p^Rp), in the system and reservoir, 
respectively. These density distributions are related to the corresponding polymer size 
distributions, P(R p ;r] n ) and P r (R p ), according to 

P P (R P ;Vn) = Pp(Vn)P(R P ;Vn) and p r p (R p ) = p r p P r (R p ) , (B.3) 
which are normalized such that 

/•oo roo 

P P (Vn) = J q dRppp(Rp]7] n ) and P r p = J Q dR p p r p (R p ) (B.4) 

are the average polymer densities in the system and reservoir, respectively. 

In free- volume theory, the polymer free energy is approximated by the sum of the 
free energy of an ideal gas of polymers confined to the free volume, i.e., the volume 
not excluded by the nanoparticles, and the entropy of a polymer coil due to internal 
(conformational) degrees of freedom: 



POO 

f p = dR p p p (R p ]7] n 

J 



ln ( P p(Rp,VnK j _ , /;(iiis( 



(B.5) 



where a(R p ; r] n ) is the free- volume fraction of polymers of radius R p amidst nanoparticles 
of volume fraction r] n and f CO nf{Rp) is the conformational entropy of a polymer coil. Now 
equality of polymer chemical potentials in the system and reservoir implies 

P P (R P ; Vn) = p r p (R P )a(R p ; r] n ) . (B.6) 

Furthermore, assuming the conformational entropy of a polymer in the system is the 
same as in the reservoir, we have 

fcorf(Rp) = - In P r (Rp) ■ (B.7) 
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(B.8) 



Combining equations (!B.3p - (lB.7p . the polymer free energy density is given by 

f P (Vn, P r p ) = Pp«eff(^n) [in - l] , 

with the effective polymer free- volume fraction a e s(r] n ) defined by equation ( TT3l) . 
Equations (IB.lj) and (IB. 81) together yield the total free energy density of equation ( TT2|) . 
Moreover, from equations (1B.3I) . (1B.4I) . and (1B.6I) . the polymer size distribution in the 
system is expressed by equation ffl~4l) . 

In practice, it is convenient to approximate the reservoir polymer size distribution 
[equation by the accurate, analytic ansatz of Eurich and Maass |75j : 



PM 



exp 



u , 9 a 
d- 

a u 



2uKq \ a uj 
where u = R 2 p , K = 0.015923, a = 0.0802, and d = 1.842 (see equation (13) in [75]). 
For the polymer free-volume fraction, we adapt the geometry-based approximation of 
Oversteegen and Roth |76j: 

a(R p ; 7] n ) = (1 - rj' n ) exp[-f3(pv p + ja p + kc p )} , (B.10) 

where r]' n = [1 — exp(— (3e))r) n represents an effective nanoparticle volume fraction 
for penetrable polymers [54J. Equation ( IB.lOj) . a generalization of scaled-particle 
theory derived from the fundamental measures formulation of density-functional 
theory [771 EBJ [79] conceptually separates thermodynamic properties of the nanoparticle 
hard-sphere fluid - - bulk pressure p, surface tension at a planar hard wall 7, and 
bending rigidity k — from geometric properties of the polymer depletant — volume 



v 



(A-k/3)RI = (Air/3)Rlq 3 , surface area 



mean curvature c 



AttR 2 p 



^■nR n q , and integrated 



p 



R„ = R n q. The hard-sphere properties are approximated by the 



Carnahan-Starling expressions [76J: 

3< 1+Vn + Vn 



n'3 



h 

/3k 



AttRI 

H 



^(2 



■Vn) 

rQ 



[1-V'n 



+ ln(l - Vn) 



Rn(l 



V'n) 



(B.11) 



[1] Pusey P N 1991 Liquids, Freezing and Glass Transition ed Hansen J-P, Levesque D, and Zinn- 

Justin J (Amsterdam: North-Holland) 
[2] Poon W C K 2002 J. Phys.: Condens. Matter 14, R859 (2002). 
[3] Tuinier R, Rieger J, and de Kruif C G 2003 Adv. Colloid Interface Sci. 103 1 
[4] Balazs A C, Emrick T, and Russell T P 2006 Science 314 1107 
[5] Asakura S and Oosawa F 1954 J. Chem. Phys. 22 1255 
[6] Vrij A 1976 Pure and Appl. Chem. 48 471 

[7] Cast A P, Russel W B and Hall C K 1986 J. Coll. Int. Sci. 109 161 

[8] Lekkerkerker H N W, Poon W C K, Pusey P N, Stroobants A, and Warren P B 1992 

Europhys. Lett. 20 559 
[9] Panagiotopoulos A Z 2000 J. Phys.: Condens. Matter 12 R25 



Crowding of Polymer Coils in Nanoparticle- Polymer Mixtures 



17 



Panagiotopoulos A Z 1987 Mol. Phys. 61 813 

Panagiotopoulos A Z, Quirke N, Stapleton M, and Tildesley D J 1988 Mol. Phys. 63 527 
Panagiotopoulos A Z 1989 Int. J. Thermophys. 4 739 

Panagiotopoulos A Z 1995 Observation, Prediction, and Simulation of Phase Transitions in 
Complex Fluids, NATO ASI Series C 460 463-501 ed Baus M, Rull L R, and Ryckaert J P 
(Dordrecht: Kluwer) 
Vink R L C and Horbach J 2004 J. Phys.: Condens. Matter 16 S3807 
Vink R L C and Horbach J 2004 J. Chem. Phys. 121 3253 
Lo Verso F, Vink R L C, Pini D, and Reatto L 2006 Phys. Rev. E 73 061407 
Louis A A, Bolhuis P G, Hansen J-P, and Mcijcr E J 2000 Phys. Rev. Lett. 85 2522 
Bolhuis P G, Louis A A, and Hansen J-P 2002 Phys. Rev. Lett. 89 128302 
Schmidt M and Denton A R 2002 Phys. Rev. E 65 21508 
Schmidt M, Denton A R, and Brader J M 2003 J. Chem. Phys. 118 1541 

Tuinier R, Aarts D G A L, Wensink H H, and Lekkerkerker H N W 2003 Phys. Chem. Chem. Phys. 
5 3707 

Rosch T W and Errington J R 2008 J. Chem. Phys. 129 164907 

Zausch J, Virnau P, Binder K, Horbach J, and Vink R L C 2009 J. Chem. Phys. 130 064906 
Zausch J, Horbach J, Virnau P, and Binder K 2011 J. Phys.: Condens. Matter 22 104120 
Sear R P and Frenkcl D 1997 Phys. Rev. E 55 1677 
Warren P B 1997 Langmuir 13 4588 

Fasolo M and Sollich P 2005 J. Phys.: Condens. Matter 17 797 

Meijer E J and Frenkcl D 1991 Phys. Rev. Lett. 67 1110; 1994 J. Chem. Phys. 100 6873; 1995 
Physica A 213 130 

Chou C-Y, Vo T T M, Panagiotopoulos A Z, and Robert M 2006, Physica A 369 275 

Moncho-Jorda A, Rotenberg B, and Louis A A 2003 J. Chem. Phys. 119 12667 

Vink R L C and Schmidt M 2005 Phys. Rev. E 71 051406 

Dijkstra M, van Roij R, Roth R, and Fortini A 2006 Phys. Rev. E 73 041404 

Schmidt M, Fortini A, Dijkstra M 2003 J. Phys.: Condens. Matter 15 S3411 

Fortini A, Schmidt M, Dijkstra M 2006 Phys. Rev. E 73 051502 

Vink R L C, Binder K, and Horbach J 2006 Phys. Rev. E 73 056118 

Vink R L C, De Virgilis A, Horbach J, and Binder K 2006 Phys. Rev. E 74 031601 

De Virgilis A, Vink R L C, Horbach J, and Binder K 2007 Europhys. Lett. 77 60002 

De Virgilis A, Vink R L C, Horbach J, and Binder K 2008 Phys. Rev. E 78 1539 

de Gennes P G 1979 Scaling Concepts in Polymer Physics (Ithaca: Cornell) 

Fixman M 1962 J. Chem. Phys. 36 306 

Flory P J and Fisk S 1966 J. Chem. Phys. 44 2243 

Flory P J 1969 Statistical Mechanics of Chain Molecules (New York: Wiley) 
Denton A R and Schmidt M 2002 J. Phys.: Condens. Matter 14 12051 
Caseri W 2000 Macromol. Rapid Commun. 21 705 

Ramakrishnan S, Fuchs M, Schweizer K S, and Zukoski C F 2002 Langmuir 18 1082 
Ramakrishnan S, Fuchs M, Schweizer K S, and Zukoski C F 2002 J. Chem. Phys. 116 2201 
Shaw S A, Chen Y L, Schweizer K S, and Zukoski C F 2003 J. Chem. Phys. 118 3350 
Kramer T, Schweins R, and Hubcr K 2005 J. Chem. Phys. 123 14903 
Kramer T, Schweins R, and Hubcr K 2005 Macromol. 38 9783 

Hcnnequin Y, Evens M, Quezada Angulo C M, and van Duijneveldt J S 2005 J. Chem. Phys. 123 
54906 

Moncho-Jorda A, Louis A A, Bolhuis P G, and Roth R 2003 J. Phys.: Condens. Matter 15 S3429 

Doxastakis M, Chen Y-L, Guzman O, and de Pablo J J 2004 J. Chem. Phys. 120 9335. 

Doxastakis M, Chen Y-L, and de Pablo J J 2005 J. Chem. Phys. 123 34901. 

Schmidt M and Fuchs M 2002 J. Chem. Phys. 117 6308 

Bolhuis P G, Meijer E J, and Louis A A 2003 Phys. Rev. Lett. 90 68304 



Crowding of Polymer Coils in Nanoparticle- Polymer Mixtures 



18 



[56] Belotelov V I, Carotcnuto G, Nicolais L, Pepe G P, and Zvezdin A K 2005 Eur. Phys. J. B 45 317 

[57] Kramer T, Schweins R, and Huber K 2005 Macromol. 38 151 

[58] Le Coeur C, Demc B, and Longeville S 2009 Phys. Rev. E 79 031910 

[59] Tuteja A, Duxbury P M, and Mackay M E 2008 Phys. Rev. Lett. 100 077801 

[60] Sen S, Xie Y, Kumar S K, Yang H, Bansal A, Ho D L, Hall L, Hooper J B, and Schweizer K S 

2007 Phys. Rev. Lett. 98 128302 

[61] Fujita H and Norisuyc T 1970 J. Chem. Phys. 52 1115 

[62] Yamakawa H 1970 Modern Theory of Polymer Solutions (New York: Harper & Row) 

[63] Eisenriegler E, Hanke A, and Dietrich S 1996 Phys. Rev. E 54 1134 

[64] Hanke A, Eisenriegler E, and Dietrich S 1999 Phys. Rev. E 59 6853 

[65] Woodward C E and Forsman J 2008 Phys. Rev. Lett. 100 98301 

[66] Forsman J and Woodward C E 2009 J. Chem. Phys. 131 044903 

[67] Vallcau J P 1998 J. Chem. Phys. 108 2962 

[68] PotofT J J and Panagiotopoulos A Z 1998 J. Chem. Phys. 109 10914 

[69] Errington J R 2003 J. Chem. Phys. 118 9915 

[70] Hynnincn A-P and Panagiotopoulos A Z 2008 Mol. Phys. 106 2039 

[71] Hansen J-P and McDonald I R 1990 Theory of Simple Liquids (New York: Academic) 

[72] Widom B J 1963 J. Chem. Phys. 39 2808 

[73] Frenkel D and Smit B 2001 Understanding Molecular Simulation (London: Academic) 

[74] Allen M P and Tildesley D J 1987 Computer Simulation of Liquids (Oxford: Oxford) 

[75] Eurich F and Maass P 2001 J. Chem. Phys. 11 7655 

[76] Oversteegen S M and Roth R 2005 J. Chem. Phys. 122 214502 

[77] Rosenfeld Y 1989 Phys. Rev. Lett. 63 980 

[78] Rosenfeld Y 1994 Phys. Rev. E 50 R3318 

[79] Schmidt M, Lowcn H, Brader J M, and Evans R 2000 Phys. Rev. Lett. 85 1934. 



